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1 Introduction 



Strongly interacting gauge-field theories, where perturbation theory is of no help, have been a long- 
standing problem. The theory of strong interactions, quantum-chromodynamics (QCD), is probably 
the most well-known theory of this type. This paper is concerned with one approach to such prob- 
lems, the Hamiltonian description of the lattice version of gauge-field theory (LGFT). For most 
investigations in this area this is not the method of choice; the majority of calculations on such 
theories are done within the Lagrangian formalism. The main reason for this is the fact that, at least 
in principle, the Lagrangian method is elegant and simple. It is based on a path-integral approach 
to the imaginary-time propagation. Much work has been done to improve the accuracy of the Monte 
Carlo method used in the numerical implementation of the Lagrangian approach, further increasing 
the viability of the method. 

Thus, anyone wishing to pursue a different path needs first to build a strong case for the merits 
of their approach. For the form of the Hamiltonian approach we shall consider in this paper we see 
four important advantages: 

1. The Lagrangian approach, being based on an imaginary time evolution, does not allow easy 
access to the vacuum wave functional. By contrast, in the Hamiltonian approach such a wave 
functional is at the core of the calculation, and we cannot avoid to calculate it. Once the 
vacuum wave functional is known, most properties of theories such as QCD, like confinement 
and chiral symmetry breaking, should follow automatically. 

2. Time-dependent phenomena can only be discussed in a real time (Hamiltonian) setting. 

3. The physical interpretation of the variables is more transparent. For example, electric and 
magnetic operators have their classical meaning. 

4. The Hamiltonian approach can studied not only with numerical, but also with analytical tech- 
niques, so we can study low- lying excitations with a harmonic approximation and we can 
disentangle the dependence of observables on the parameters in the approximation. 

These advantages are accompanied by important disadvantages. Due to the fact that the Hamil- 
tonian requires a quantisation plane to be defined, we lose the manifest Lorentz symmetry, that arises 
so naturally in the Lagrangian approach. In the full solution to the Hamiltonian problem this should 
be recovered, but approximations can introduce spuriosities. Furthermore we have to fix the gauge, 
and the degrees of freedom in the Hamiltonian formulation of the theory will then naturally depend 
on the choice of gauge. The issue of gauge fixing is a rather non-trivial exercise, as we shall see in 
this paper. 

We shall first, in Sec. ^, give a short introduction to our version of lattice gauge theory. In 
Sec. |3| we then discuss how we can use the idea of a maximal tree to define a fully gauge-fixed 
Hamiltonian, which is derived in an explicit form. In Sec. ^ we discuss those constraints relating 
the trace of powers of SU{N) matrices, and we define a symbolic method to construct eigenvalues 
of the electric energy operator, which should be useful in any approach to Hamiltonian LGFT. We 
then study the one plaquette problem in Sec. ^, first in its generality, showing how we can derive 
a simple A^-fermion problem equivalent to the SU{N) one-plaquette problem. We discuss both the 
weak-coupling (harmonic) and the large-A^ limits. Explicit results for the one plaquette spectra as 
function of the coupling constant are shown for U{1), and SU{N) LGFT's, with 2 < < 5. Finally 
we discuss the results and future directions in Sec. |^ 

2 Hamiltonian lattice gauge theory 

The degrees of freedom in a gauge theory, such as electromagnetism, can be divided into two sets. The 
first set consists of the charges and currents in the system together with the associated static electric 
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and magnetic fields. The second set contains the electromagnetic waves. However, since a charge 
in a moving frame is a current, relativistic coordinate transformations mix electric and magnetic 
variables. Moreover, since the velocity of electromagnetic waves is constant, the two polarisation 
directions, perpendicular to the propagation direction, also undergo complicated transformations 
under a change of frame. Therefore a relativistically invariant formulation of electromagnetism is 
very useful. 

As is well known, the Lagrangian of electromagnetism, without static sources. 



C 



\j:{Ei-B^) , (1) 



can be formulated in terms of a covariant four-vector potential as 

C = ^{d^A,-d^A,)^ , (2) 

where Bi = eijkdjAk and Ei = diAo — doAi. We use a notation in which roman indices i,j, k, . . . 
run over the spatial, while greek indices . . . run over the values 0, 1, 2, . . . ,D, where D is the 
number of space dimensions and the temporal axis has index 0. The introduction of the vector 
potential, or gauge field, serves a number of purposes. Firstly, since the gauge field transforms as 
a vector, the Lagrangian is now manifestly relativistically invariant. Secondly, it is formulated in 
terms of canonical variables A^ and d^A^ which should make quantisation more straightforward; 
and thirdly, the two constraints for the electric and the magnetic fields, the source-free Maxwell 
equations diBi = and doBi = —eijkdjEk, are automatically satisfied. The price we pay for this 
is the appearance of superfluous degrees of freedom, the gauge, and an entanglement of the sources 
giving rise to these variables. 

The original Abelian gauge theory of electromagnetism was extended by Yang and Mills |j^] using 
gauge fields of more complicated structure, including internal degrees of freedom. This generates 
self-interactions since the gauge fields do not commute, but are chosen to obey the commutation 
relations of a specific Lie algebra. We shall concentrate on a gauge-field Lagrangian where the field 

is an element of the Lie algebra su{N), 

£ = ^Tr[F,,Fn . (3) 

Here the field tensor is the skew-Hermitian matrix i^^^y = d^Ai, — d^A^ — i[A^,Ay\, and the field 
variable A^ is an element of the algebra, conveniently parametrised as 



A.^g-X'^Al . (4) 



The A"^ are the A^^ — 1 generators of the Lie algebra, satisfying the commutation relations 

[A^ A''] = 2ir''^A^ . (5) 

The index a thus runs from 

1 to A^2 _ i_ The A"" can be represented by traceless N x N matrices, 
normalised such that their squares have trace 2, as can be seen from the anticommutation relations 

{X-,X'} = 2d'^'''X' + ^6abI . (6) 

In Eq. (^) we have absorbed the coupling constant g in the field A^, so that we can interpret the fields 
geometrically, since the field tensor is now the curvature that follows from the covariant derivative 

d^^d^-i[A^,.] , (7) 
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I.e., 



I-IVI 



(8) 



As we are interested in the Hamiltonian, we perform the standard equal-time quantisation and 
reformulate the Lagrangian in generalised electric and magnetic fields. This is strictly speaking a 
3-dimensional result, since this interpretation requires the use of three-dimensional algebra. We shall 
nonetheless use the result below for other dimensions as well. We find 



C 



1 



Tr 



D 



.fc=i 



■Bl) 



(9) 



where Bi = —^^^ijkPjk and Ei = FiQ. Since we wish to impose the temporal gauge Aq 
separate the Hamiltonian in two parts, thereby isolating the Aq dependent part: B 



0, we 



£ = ^Tr 



D 



Bl] 



.k=l 



+ ^TT[Aog + AoXiAo)] 



(10) 



where we added a total divergence. The function X is second order in ^Iq and does not contribute 
to the equations of motion, or the constraint equations in the temporal gauge {Aq = 0) discussed 
below. Since the time-derivative of Aq does not occur in the Lagrangian, the (Lagrangian) extremal 
action variational principle leads to an equation of motion for Aq that is a time-independent algebraic 
equation, which shows that the Aq takes on a time-independent constant value. This set of equations 
(one for each colour index) are the non-Abelian analogue of the Gauss' law constraint, which in the 
absence of colour charges take the simple form 



a'^(x) = 



ill) 



where 



D 



D 



g'^(x) = ^[d,Et{^)+gr'^AU^)Efi^)] = ^d,Et 



(12) 



i=l 



i=l 



The components of the fields can be obtained via the relation 



(13) 



The constraints obey the same commutation relations as the generators of the gauge group. Thus 
Gauss' law cannot be implemented as a strict operator condition as it leads to contradictions, since the 
non-commuting constraints cannot all be diagonalised simultaneously. However, within the physical 
(in this case colourless) subspace defined by 



g"(x)|Phys) = 



(14) 



no such problem arises, since the eigenvalue of the commutators is also 0. The space of states consists 
of wave functionals, which are functionals taking values on the SU{N) group manifold. We find A^'^ — 1 
functional conditions on each wave functional, consisting of functions on the group manifold at each 
space point. 

As is well known, quantisation of problems involving redundant degrees of freedom (i.e., where 
some of the equations of motion are constraints) is quite involved. The two main techniques used 
are Dirac and BRS quantisation, and they require a large amount of additional analysis. For more 
details one can consult the seminal work by Dirac ||3|,|3, as well as Refs. [0,|5|- 
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If we are able to work within the physical subspace only, one can ignore these formal problems 
and define the quantisation of the canonical momenta 11^ = (9o^a by 



which involves a functional derivative ||, 10| with respect to the field variables. 

Since Aq is not dynamical, we cannot associate a canonical momentum with it. We therefore use 
the temporal gauge, Aq = 0, which leaves us with residual gauge freedom 0(x) independent of the 
time coordinate, i.e., under the transformation 

^^^(x) = 0(x)A^(x)</.-i(x)+^[5^</.(x)]0-Hx) , (16) 
^F^,(x) = </.(x)F^,(x)<^-i(x) , (17) 

where ^ G SU{N), the Lagrangian is invariant. 
2.1 Discretisation 

Field theories suffer from singularities, both in the infrared and ultraviolet limits. In many interesting 



cases, such as QCD |11|, these are renormalisable. Rather than dealing directly with the continuum, 
we shall regularise the problem by introducing a simple hypercubic lattice in the D-dimensional space, 
with lattice spacing a. Since we are pursuing a Hamiltonian approach, time will remain continuous. 
In this paper we shall concentrate on pure gauge theory, without explicit charges (i.e., quarks). In 
this case the systems is described by a set of gauge fields Af at each point of the D-dimensional 
lattice. It looks like these could in principle carry D{N'^ — 1) degrees of freedom, as the group 
SU{N) has N"^ — 1 generators A". However, even after restricting the gauge fields by imposing the 
temporal gauge Aq = 0, the Lagrangian is still invariant under a limited set of gauge transformations 
that do not violate the temporal constraint. Thus the apparent number of degrees of freedom is still 
larger than that required to specify the dynamics uniquely. These superfluous degrees of freedom 
arise since we have not fully exploited the Gauss' law constraints. 

In order to analyse this problem fully we first investigate the structure of the Hamiltonian on the 
lattice. The gauge fields 

= ^2 E ' (18) 

a=l 

are Hermitian, since A" is Hermitian. In the Schrodinger (wave function) representation the effect 
of these fields can be incorporated as a non-Abelian change of phase of the wave function between 
different points, a simple generalisation of the Aharonov-Bohm effect in QED. This geometric inter- 
pretation can be represented by a group element, a unitary transformation summing all the small 
phase changes along the path connecting the two points, 

U{^i,^2) = Pexp\i [ ' A^{^)dxA , (19) 



where we have introduced a path-ordered product, denoted by P, since the quantities .A^(x) for 
different values of x are non-commuting matrices. 

In practice we will only use this unitary transformation (which for obvious reasons is also called a 
parallel transporter) on a link between two nearest-neighbour lattice points i and j where = xj — xj 
is a primitive lattice vector in the direction /i (we shall use the bracket notation () to denote nearest- 
neighbour pairs) 

C/(ij) = Pexpl^ij A^(x)d3;^| = exp {za^^^(x(ij))} , (20) 
C^Cii) = exp{-m/,A^(x(ij))} = C/Jj^ . (21) 
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(0,a„) 




(0,0) (^x'O) 
Figure 1: The path for the elementary plaquette. 



In Eq. (2C) we have defined at the midpoints of the hnk, 



(22) 



Since is an element of the su{N) Lie algebra, [/ is a SU {N) matrix, and satisfies U^^ = and 
detU = 1. The quantity j4^(x^ij^) is a geometrical average of along the link (ij) in direction 
Since all future discussions will only concern this averaged field, we suppress the bar from here 
onwards. 

The four links around a primitive square (usually called plaquette) on the lattice will give, when 
we take the trace, the simplest gauge-invariant quantity on the lattice, and for that reason is often 
used in lattice problems. As an example we take a plaquette in the xy plane, (0,0) (ox^O) 
{ax,ay) (0,ay) —>■ (0,0), see Fig. |l|. This leads to (the notation aj is used as a shorthand for 
times a unit vector in the direction i) 



Uo 



exp {iaxtty {AyAx{m) - Ar,Ay{m) - [Ax{m),Ay{m)]) + 0{a^} 



(23) 



where m = (a^^ + ay)/2 is the centre of the plaquette, and we have assumed that both and Uy are 
of order a. We have used the Campbell-Baker-Hausdorff formula to combine the exponentials. 



exp(^) eMB) = e^p!^A + B + ^[A,B] + -^[A,[A,B]] + ^[[A,B],B] + 

The operator Aj is the lattice derivative, defined as 

Aj{m + ai/2) - Aj{m - a,i/2) 



AiAj{m) 



Eventually we find, using 



U 



-a'^iF^y+0(a^) 



(24) 



(25) 



(26) 



We would like to associate the trace of U with a lattice version of the QCD Lagrangian. The 
requirement is that in the a — > limit the Lagrangian becomes the continuum one. Thus we must 
analyze the expansion of the trace. If we expand U\j in a Taylor series in a, we find that the and 
terms are traceless, and so are all the terms of the order except for the term —i{Fxy)'^- Thus 
it seems natural to define TrC/^ as the lattice potential. Unfortunately, this is a complex quantity, 
and we need to combine it with its complex conjugate to get a real result. 



1 



1 



■Ti[2-Uu-U^^] + 0{a 



(27) 



Clearly this approaches the continuum result as a — > 0. This is the one-plaquette contribution to the 
magnetic part of the Wilson action [p^ ,|l3|. 
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In a Hamiltonian formalism we need to construct canonical momenta. The natural choice is the 
differential realisation of the generators at each lattice point (i.e., the generalised angular momentum 
operators for the group SU{N)). These act naturally on the original lattice variables A^{xi), but 
as stated above we really would like to formulate the problem in terms of link variables C/^ij^ . Using 
the definition of f/^ij^ as a path-ordered product, Eq. (|l9|), we find that the action of the momenta, 
conventionally called electric fields E, is given by 

£;^(xi)C/(ji) = ic/(ji)A» , (28) 

where j = i + e^, so fj, must point along the link (ij) in the positive direction. The second relation 
follows from the Hermitian conjugation. This can also be written as a differential operator acting on 
the link-averaged gauge fields AJJ(x^ij')), 



We shall also use the notation Ef^ for this same operator. 

Ignoring problems with overcompleteness, we can derive the Kogut-Susskind Hamiltonian, 



(ij) 

where n is the length of the lattice, and D is the number of space dimensions. The quantity U (rrii^fc^/) 
is the Wilson plaquette around the face centre rrii in the directions k and I, 

U{uii^k,l) = t^i,(i+afc)^(i+afc),(i+afc+aO^(i+afc+aO,(:+a!)^(i+a;),i ' (^1) 

where mj = Xj -|- ^(a/j + a;). We have written down the slightly more general form of Hamiltonian 
obtained when we realise that the coupling constants gt and gs for the spatial and time components 
of A can be different, since the Hamiltonian breaks manifest Lorentz invariance. In principle one has 
to solve a renormalisation problem to find the physical couplings that correspond to these bare ones, 
a problem that we shall not address in this paper. 

2.2 Gauss' law 

The constraints, Eqs. (|ll],||), can also be interpreted as arising from the invariance of the gauge-fixed 
Lagrangian under residual gauge transformations <j>. Unfortunately, the result (^) is only correct for 
continuum theories, since it involves the field rather than its lattice analogue U. The derivation 



of Gauss' law constraints on the lattice is discussed in Ref. |14]. The approach is based on the fact 



that a wave functional, which dependens on the overcomplete set of link variables, 

^{{U}) , (32) 
Should be invariant under a gauge transformation on a single lattice site, 

- e^'^'C/i,i+e,e'^'+^M (33) 
C/i_. ,i - e^^-^^Ui.e,,i(^'^ . (34) 
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Table 1: The number of degrees of freedom on a hypercubic lattice of size n x n x ■ ■ - n, for different 
numbers of space dimensions. The D = 1 case comprises n plaquettes on a line 



dimensionality 




2 


3 


D 




number of sites 


2(n+l) 


(n + l)2 


(n+ 1)3 


(n + i; 




number of links 


3n + l 


2n(n + 1) 


3n(n+ 1)2 


Dn[n - 


f 1)^-1 


number of plaquettes 


n 




3n^ 


\D{D 


-l)n^ 



The invar iance can be written in terms of generators as 

D 



exp 



where are the left-handed generator for SU{N), 



*({[/}) =0 



-A" 



(35) 



(36) 



which give rise to the right-multiplication in Eq. (^). 

Prom the infinitesimal from of Eq. (35 we find the lattice form of the Gauss' law, 



D 



(37) 



It is easy to show that 



p{L)a 



2 



1,^ ' 1,11 J 

We thus find the lattice version of the Gauss' law generators 



(38) 
(39) 
(40) 

(41) 



In the limit that the lattice spacing vanishes this simplifies to the original Gauss' law generator, 
Eq. (0). 

From Eq. (|4l] ) we see that we can define a covariant lattice derivative, which takes the matrix 
form 

■Dij,/. = 5i+e^,j - Dab{Ui^e^;,)Si^e^^} ■ (42) 

This construction of the Gauss' law makes it clear why it is so attractive to work with functions of 
traces of the product of U along closed loops (Wilson loops), since those varaibles are automatically 
gauge invariant. 

Thus the number of constraints equals the number of group generators times the number of sites, 
and the number of degrees of freedom equals the number of generators times the number of links. 
The number of unconstrained degrees of freedom is clearly the latter number minus the former. As 
shown in Table |l], for one and two space dimensions, in the limit of infinite extent, this is exactly 
the number of plaquettes times the number of gauge degrees of freedom. This is not true in higher 
dimensions, and we shall have to do some more work to get around this problem. 

Another interpretation of this procedure is that Gauss' law fixes the longitudinal electric field. 
This field is a direct consequence of the charge density distribution in the system. Since we will work 
in the colourless sector, the longitudinal field is zero, in agreement with our results above. 
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3 Explicit Hamiltonian 



The link variables are a sufficient set of variables, and, when combined in closed contours, they are 
gauge invariant. This means that projections on the colour-neutral sector by integration over the 
gauge group, as seen in many Lagrangian lattice gauge approaches to the problem, have no effect. 
The price we have paid for the Hamiltonian approach is the breaking of explicit Lorentz invariance, 
and the problem of determining the physical subspace. It is also easy to see that the collection of 
all possible contour variables form an overcomplete set. It can be shown explicitly that there are 
relations between different combinations of contours, and also between contours winding a different 
number of times around the same path. 

For a systematic Hamiltonian approach, we want to define an exactly complete (i.e., neither 
overcomplete nor incomplete) set of variables, and we must know the effect of the electric and 
magnetic operators upon these variables. Part of the overcompleteness is due to the fact that we 
still have residual gauge degrees of freedom. Therefore we shall fix of all attempt to fix the gauge 
fully. Some natural choices, such as Coulomb and maximal Abelian gauges |14-16|, are not suited to a 
Hamiltonian approach. The Coulomb gauge is highly non-local, which makes localised approximation 
hard to construct, and the maximal- Abelian gauge requires a diagonalisation procedure for each value 
of the field, not an attractive procedure in an operatorial approach. The Coulomb gauge was applied 
in the Hamiltonian framework by Christ and Lee [17|, but that approach has not shed any light on 
the properties of QCD. 



Inspired by initial work of Miiller and Riihl [18|, and Bronzan |ll9|-p^], we shall use a maximal-tree 
gauge where we gauge-fix a maximal number of links, which leaves a complete set of variables in the 
Hamiltonian. We combine these link variables to form closed contours, which is useful in the absence 
of colour charges. This is still not enough, since another form of overcompleteness creeps in related 
to the number of degrees of freedom in traces of SU{N) matrices. In the next section. Sec. |^, we 
investigate the question of what conditions exist among different trace variables. 



3.1 Conditions on variables 

As in electromagnetism the magnetic field should be divergenceless. This is a direct consequence of 
the skew-symmetry of the field tensor, 

D 

^A(m)S,(m) =0 , (43) 

i 

which follows from the Jacobi identity concerning the vanishing of the sum over all cyclic permutations 
of the double commutator X^cyciicl-^i' i^j^^k]] = 0. Therefore the magnetic fiux through a closed 
surface should be zero. Since we have contours and more than one component of the magnetic field 
in non-Abelian SU{N) gauge theory, similar relations can be derived for arbitrary contours a. In 
order to combine the magnetic fiux through a set of contours we define |2j] 

< = Ti-m , (44) 
= TriX'^U^] . (45) 

The set {VKq,!^^} which uniquely determines Ua, as can be seen from the decomposition Ua = 

From commuting the matrices in a trace of arbitrary length we can derive relations between 
contours and their decomposition in smaller contours. For example for SU{2) where is real, and 
is purely imaginary we find 

Txp^Up] + TV[[/„C/. 1] = Tr[C/,] Ti[Up\. (46) 
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For a general SU{N) we find for the trace of three matrices that 



2 



Tr[UaUpUy]+TiiUaU^Up] = 2^—^TT[Ua]Tv[Uf,]Tv[U^ 



Ar2 



+^ ^ d'^^'^VF^Ty^W^^ . (48) 

a,b,c 

For the case of SU{2) and SO(3) these equations simplify, since in that case d"'^'^ = and f°'^'^ = e°'^'^. 
These relations define a set of constraints, usually called "Mandelstam constraints" |24|. 

Not only do relations exists among contours, but also for paths winding around the same contour 
more often than — 1 times. These trace along such a contour can be re-expressed in terms of the 
first — 1 trace variables, 

Tr[C/"]; n= 1,2,--- ,Af- 1 . (49) 

So we find 

Tr[[/P] = X] Tr[C/]«i Ti[U^f^ ■ ■ ■ T^[U^-^f^-^ , (50) 



where ^j=i ^ P- Generally we do not need to know these relations for all values of p. If we 
determine the eigenstates of the electric operator the necessary relations follow. We will derive the 
values for for some useful cases in Appendix A. 

For SU{2), which has only one independent trace variable Tr[C/], the relations are, as can be 
derived from the representation, 

Tt[U-^] = 2cos(j) = Tr[U] , (51) 



and the lowest-order relations for SU (3) are 

Tr[U-'] = lTr[[/f -iTV[C/2] , (53) 
Tr[[/3] = 3 + ^Tr[C/2]Tr[C/] - ^Tr[C/]3 . (54) 

3.2 The gauge-fixed Hamiltonian 

From the counting of degrees of freedom in Table || we know that the set of all link variables must in 
general be overcomplete. Since the link variables are still one of the most attractive sets to use, we 
choose to fix the gauge as much as possible, based on some ideas developed by Miiller and Riihl [18|, 
and Bronzan and collaborators [|l9|-p3|. The approach is based on fixing the gauge along the links 
of a maximal tree, a collection of links on the lattice which connect any two lattice sites by one and 
only path. 

{Uij} ^ {V^i,W^^\W e tTee,V ^ tree} . (55) 

These W variables can serve as our "longitudinal electric fields", and are determined by the 
gauge fixing. The links not on the tree form the base for the "magnetic variables". However, as 
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Figure 2: A maximal tree as defined in Eq. (|55|). 

they stand, these variables are not invariant under local gauge transformations, which means that 
the wavefunction cannot depend directly on these variables, because that would contradict gauge 
invariance. Hence we combine them with a path over the maximal tree from and to the origin, 

Xii = Wo,iViiW^j , (56) 

see Fig. ^. In order to have a unique way of defining the links on the maximal tree, we have oriented 
the links on the maximal tree such that they point along the direction away from the origin, which 
is why an inverse appears on the far right in Eq. (^), since Wqi is defined as the product along the 
maximal tree, 

7(i)-l 

Woi= n ^Ju.+i > (57) 

k=0 

where 7(i) is the number of links on the maximal tree between the origin and site i. Note that there 
is no trace in Eq. (|56|), and thus all the variables transform in the same way under local gauge 
transformations with the gauge transformation at the origin, and are invariant under all other local 
gauge changes. We know that when we fix the gauge we cannot fix a global gauge transformation, 
and we are thus led to identify this with the one at the origin. 

We can both use and Xji as variables, but since these are clearly related by an inversion, we 
shall just use one of them. In the sum over links in the electric operator we shall also just use one 
order, which we choose again to be such that the link point along a branch of the maximal tree away 
from the origin for tree links, and such that the first non-zero component of the direction of the link 
is positive for the remaining ones. 

The electric operators E?;y defined on every link naturally fall into two classes: those acting on 
the maximal tree variables W and those acting on the variables V. Let us first look at those acting 
on a tree variable, and evaluate the action on a variable Xij. Again we distinguish two cases, since 
the link can be contained in either (actually even in both!) the path from to i, or in the path from 
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Figure 3: A sketch of the path in a typical variable X. 



j to 0. If the electric operator acts in the first part we find 

7(i')-l 7(i)-l 

^i'j'^ij = n (^i'j'^i'j') n ^j.j.+i^ij^oj' 

fc=0 7(j') 



= -^W^oi'A^T^oT'^iJ • (58) 
Similarly, being much less explicit, when acting on the last part of X we have 

^i'j'^iJ = ^^ijW^oT'AWoi' • (59) 

Finally the action on the non-tree variable will also lead to left-multiplication by A", due to our 
choice of ordering, and we have 

Aj^iJ = -^^oiA^W^oT'^iJ ■ (60) 

These expressions are clearly highly suggestive of defining a new operator that leads to left or right 
multiplication by ^WoiA^Wj^^. Using some simple elements of representation theory, we actually 
know that 

W-^X^W = J2D"^iW)X'' , (61) 
b 

with 

D^\W) = ^Tv[W~'^X''WX''] . (62) 
If we define a set of "intrinsic" electric operators by 

= L>«^(Woi)4 , (63) 
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it is a matter of straightforward algebra to show that these new operators act much more simply 
on X, leading to a left multiplication by — if the labels on the link referred to by the electric 
operator occur in the same order in the product, and a right multiplication by ^A" if it appears in 
inverse order. Now we make use of the unitarity of the representation functions, 

Af2-1 

DlbDac = he , (64) 

a=l 

to show that the square of the electric operator, which is the expression that occurs in the Hamilto- 
nian, is invariant under the transformation, 

a a 

If we assume that the wave function can be expressed in terms of the variables X^^ , we need to 
find the action of the sum over the whole lattice of the electric energy, proportional to the square 
of the electric field, on this function. Since there are two electric operators in the square, they can 
either both act on the same variable, or on different variables. In the first case we find contributions 
where is multiplied from the left with — ^A", for all electric operators along the path leading up 
to the site i, as well as the link (ij) itself, and multiplied from the right with ^A" back from the site 
j to the origin, 

{£^,yYx,i = -i(7(i) + l)X'X,i + ^7(j)^ijA" . (66) 

(i'j'> 

In the second case each electric operator acts on a different variable X, 

Y [^i'j'^ij] [^ri'^ki] = I (7(i nl) + <5(ij),(ki)) A-^XijA^Xi™ 

(i'j') 

-^7(inm)A«XijXi™A« 
-i7(jnl)XijA"A'^Xi„ 

+^7(j n m)XijA«Xi„A'^ , (67) 

where the notation 7(i n j) denotes the number of links common between the maximal-tree paths 
from to i and from to j. 

We now wish to show in what sense the variables X are sufficient. Given a general wave function 
in terms of both the X and the W variables we investigate when this wave function lies in the 
colour-neutral sector. Then it has to satisfy Gauss' law, 

g-{i)\^) = 0; Vi^O , (68) 

which can be reformulated as a local gauge transformation generated by Gauss' law, 

e^^^e-^WlV;) = l^'^^V) = IV') . (69) 

In a {X, W} representation the wave function changes to 

{{x^j,wi^}\U) = {W'Ho)x^,im,(i>~\min.H^m) ■ (70) 

We have already identified the transformation <^(0) as the global gauge transformation, and we are 
not concerned about non-invariance under these transformations. The W variables are clearly not 
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invariant under local gauge transformations, and a function in the physical subspace is thus seen to 
be a function of the X variables only. 

This approach can also be identified with a gauge fixing where we set the matrices on the maximal 
tree links to unity, 

{Vij,VFki}-{Xij,l} , (71) 

but we keep track of the fact that the electric fields on the gauge-fixed links still act non-trivially on 
the physical variables. 

Let us now study the effect of the global gauge transformation 0(0), 

|^WV)=f/(<^)lV') • (72) 

We know that a colourless wave function transforms as a singlet under global gauge transformations 
(i.e., it is invariant under <^). One way to make sure that the wave function satisfies this condition 
is to combine the X variables in the wave function in traces, each of which transforms as a colour 
singlet, 

TrK], TrKX|], Tr^X^X^], (73) 

However, as we have seen before, this set is overcomplete; there are constraints among the different 
trace variables. In another approach to this problem we can relate the variables to group elements 
in colour space, and in the colourless sector the wave function can depend on all scalars that can be 
defined with these group elements. If we use the Mandelstam variables [^4| 

H^„ = Tr[X„]; W^ = Tt[\''X^] , (74) 

we find that there must be a constraint among them, as there are only A^^ — 1 independent degrees of 
freedom. Any unitary matrix V = (vi, V2, • • • , vat) G SU{N) has 2N'^ degrees of freedom satisfying 
A^^ + 1 constraints {deiV = 1, Vj • Vj = 5ij). Although the number of degrees of freedom equals 
the number VF^'s, the constraint do not fully fix the additional degree of freedom [Wa), which can 
take on several discrete values, corresponding to discrete gauge transformations. For example, SU{2) 
can be represented by a unit vector on the four-dimensional sphere. In this case the discrete gauge 
transformation is the sign of W^, which determines whether the vector is on the positive or the 
negative hemisphere. We shall study an approach to these problems in the next section. 

We shall close this section by making a particular choice of maximal tree, and defining the 
Hamiltonian within the physical subspace explicitly, in a recursive way. First of all (for D = 3) we 
choose all links on the x-axis for y = z = 0, all the links in the y-direction for z = 0, and all links in 
the z direction, as in Fig. ||. In a more mathematical notation the maximal tree can be defined as 
the union 

[J{[iix,0,0){ix + 1,0,0)]; [iix,iy,0){ix,iy + 1,0)]; [{ix,iy,iz){ix,iy,iz + 1)]} • (75) 
The overlap of two paths is given by 

7(1 nj) = mm{ix,jx) + Kj. min(ij^,iy) + ^i^jji^j^ Ta.\n{i^,j^) . (76) 
We start from the {D = 3) Kogut-Susskind Hamiltonian, Eq. (|30|), and let it act on a function of 
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the X variables. For the electric part we have 



He 



(kl),(ij>0ree a 

-7(1 ni) 
-7(knj) 

+7(1 nj) 



A"X 



JXki 

5 



A"X 



(JXki 

5 



5Xki 

5 ' 
(5Xki 



XvX 



u J 



A"X 



u J 



(77) 



where the functional derivative combined with the matrix A" is meant to denote the insertion of A" 
at a position where the variable occurs, e.g., 



A«X-, 



5 



6Xi 



fixf,.,) = 6^,6^, (fix^x^) + f'ix.^x-xfr') + ... + f\xfr'x-x,^) 



(78) 



In the magnetic part we can set all of the links on the maximal tree to unity, 



Hm = J2 ^ 

plaquette 

= EH 



2 i^plaquette ^plaquette 



!^tree— >1 



2-X, 



(i,l,0)(i+l,l,0) ^(i,l,0)(i+l,l,0) 



2-Xf. 



(jJ,0)(i+lj,0)-^(ij+l,0)(i+lj+l,0) ^{jJ+l,0){j+lj+l,0)^(ij,o)(i+l,i,0) 



+ X] '^^ ~ ^(iJ>fe)(j+lJ,fc)^(i,j,fe+l)(i+lj,jt+l) ~ -^{iJ,k+l){i+l,j,k+l)X(ij,k){i+l,j,k) 



+ 



2 - Xr, 



(i,o,i)(i+iAi) ^(i,o,i)(j+iAi) 



2-X, 



i 

+ E Tr 

i,j,k 

(j,fe)7^(0,l) 

+ ^ TV [2 - - X-i. ^ 



-1 



X, 



^x; 



+ 



i,3,k 



2 



-X, 



(iJ,k)(iJ+l,k)X{iJ+l,k){i+lJ+l,k)Xl^illj^k){i+lJ+l,k)^{iJ,^^ 



j,k) 



(79) 



The first sum corresponds to plaquettes is the xy-plane, touching the x-axis, which contain only one 
link not on the maximal tree. The second sum contain all the other links in the xy-plane. The third 
and fourth sum contain plaquettes in constant-y planes, and the fifth and sixth in the constant-x 
planes, respectively. The last sum contains the plaquettes in the z-planes for z 0. 

3.3 Boundaries 

Let us spend a little time counting degrees of freedom, to see how close we are to the expected number. 
On a finite lattice of length n in D dimensions there are (n + 1)^ lattice points, and D{n + l)^~^n 
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Table 2: The number of degrees of freedom in the maximal tree approach lattice of size nx nx • • • n, 
for different number of space dimensions. 



dimensionality 


2 


3 




D 


number of sites 


(n + l)2 


(n + 


1)3 


(«. + 1)^ 


number of links 


2n(n + 1) 


3n(n 


+ lf 


Dn(n + 1)^-1 


number of links on tree 


n{n + 1) + n 


3(nH 


-1)2-1 


(n + l)^ - 1 


number of independent variables 




2n3 - 


f3n2 


{D - l)(n + 1)0 - D{n + 1)^-1 + 1 



links. The maximal tree contains n + n(n + 1) + n(n + 1)'^ + . . . + n{n + l)^"-*^?! = (n + 1)^ — 1 
links, as can easily be seen from (the Z)-dimensional generalisation of) our explicit choice of maximal 

tree. Therefore there are {D - l)N^ - N^-^ N remaining links, and a similar number of 

variables X. Each of these variables has N'^ — 1 degrees of freedom. The corresponding total number 
of freedom is thus just N'^ — 1 bigger than the number of unconstrained degrees of freedom discussed 
in Sec. |2.2| . The additional degrees of freedom are associated with the global gauge degree of freedom, 
represented by the gauge freedom at the origin. 

Our method shows some similarity to attempts to quantise QCD on a torus |^5|. On a torus we 
have to impose periodic boundary conditions, where identify the n + 1-st point on the lattice with 
the first point. The same maximal tree applies, however, we cannot fix the fields on the links from 
the nth point to the n + 1st point, since in that case we would have a closed path, and would violate 
the uniqueness of paths on the maximal tree. Thus we have an additional D{n + 1)^~^ free links 
connecting the boundaries of the lattice. Therefore there are D{n+l)^~^ extra independent variables 
X. D of these new degrees of freedom are Polyakov loops, which wind around the lattice ||26|| . 

One of the goals of our work is to apply many-body methods to the Hamiltonian, in all likelihood 
eventually we will apply the coupled cluster method, which expresses the wave function in terms 
of correlations between different variables. For an efficient implementation, we assume translational 
invariance, which restricts the number of independent coefficients in the wave functional considerably. 
Correlations on a periodic lattice are truly translationally invariant. However, on a finite lattice, if we 
take the lattice size to infinity the system should have a inner region where translational invariance 
is observed. In practice we assume an infinite lattice. 



4 Further constraints 

Even though we have drastically reduced the number of degrees of freedom, additional complications 
arise when we impose colour neutrality on the wave function. Here, the natural choice is to work with 
traces of products of the variables X, as discussed above. A suitable approach would be to construct 
a basis of eigenstates of the electric part of the Hamiltonian and calculate matrix elements of the 
magnetic energy between these states. Such an approach is a quite natural calculational scheme for 
the Hamiltonian approach. One can also use the method inherent in the Lagrangian calculations. 



based on invariant integration over the full group |27], using Monte Carlo integration. However, for 



a proper Hamiltonian approach this discards many of the advantages of the method. 

To find eigenstates of the electric operator, one can resort to three general approaches. Firstly, 
group theory gives us, in principle, a way to construct general eigenstates, the group characters. 



However, for a large basis, and N > 2, this is hopelessly involved unless it can be automated, 
and we see no easy way to do this. 

A second approach is based on integrating configurations, and constructing orthogonal combina- 
tions from them. In this case one must start off with much larger overcomplete sets of configurations 
and at increasing orders the integration, generally based on Creutz's integration method [p^, Bo|, 



tends to be more and more involved |31]. 



The third approach is based on the action of the electric operator itself, which leads to a block- 
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diagonal matrix which has to be diagonaUsed to recover the eigenstates. In combination with a 
symbolic method, this seems to be the most powerful approach, which allows one to tackle any 
arbitrary SU (N) group. We shall develop this method here. The method also applies to schemes 
other than the maximal tree method, and we shall therefore discuss it in a slightly more general 



framework (see also Refs. [32, 3^). 

We shall look at states in the colourless sector, where all variables can be formulated in terms of 
traces. Clearly only (products of) traces defined over the same set of links can mix under the action of 
the electric operator (since states with differing links are trivially orthogonal). Within such a choice 
the full set of eigenstates separates into families which are characterised by the number of lines along 
each link. In the case of a sing le plaquette we have the collections: {!}, {Tr[C/]}, {Tr[U]^ ,Tv[U^]}, 
{Tr[[/'^], Tr[[/^] Tr[C/], Tr[[/]'^}, • • • . Within each collection, we can form linear combinations of the 
elements that are eigenstates of the electric operator. However, since for SU{N) the contour Tr[?7P] 
with the winding number p > N can be expressed in terms of contours with lower winding numbers, 
one finds, as one re-expresses these contours with higher winding numbers, that some eigenstates do 
not contain leading order terms, but reduce to states of other, lower-order, collections. 

The symbolic approach starts by writing explicit indices on all the traces, labelling the column 
indices ii,... ,in- Since every member function has no free indices, the set of row indices must 
be equal to the set of column indices, but the order in which they appear can clearly be different. 
Actually, all the other members in the collection correspond to a distinct permutation, but not every 
permutation generates a new member of the collection. For example a segment with three lines has 
six permutations: 

(123), (213), (132), (321), (231), (312) , (80) 

where, in the case of a single loop, when one end of the segment is contracted with the other, the first 
element corresponds to Tr[C/]^, the next three elements to Tr[C/^]Tr[C/], and the last two elements to 
Tr[C/^]. This follows easily if one considers the general form: 

{Pl,P2,P3)= ^ Ui^ipJJi^ip^Ui^ip^ , (81) 

where P = {Pi, P2, P3) is a permutation of (ii,i2i^3)- However, in the case that two lines are con- 
nected in one loop W and the last line is connected to a separate single-line loop V, i.e. ^(W?7)(VFC/)(y[/), 
all the permutations are distinct. 

If we wish to evaluate the action of the electric part of the Hamiltonian on such a segment k), 
we first need to understand what happens if we look at a case where each electric operator acts on 
a different matrix U, 

a 

In the last step we have just changed the order in which we write the U matrix elements, suggesting 
that we can interpret this part of the action of the electric fields as the interchange of Pi and P2. 
The result Eq. ( |8^ ) is based on a standard relation for the SU{N) generators A°, 

a=l 

Using the relation (|8^ ) we find that the only non-trivial effect of the electric part of the Hamiltonian 
is an exchange of pairs, 

Y,EHiJ,k) = ^[{j,z,k) + {z,k,j) + {k,j,i)] - A fH+n{N' - 1)) {^,J,k) , (84) 

links V V / / 
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where L is the length of the segment (i.e., the number of electric operators acting on it) and n = 3 
is the number of lines. Note that the electric operator for the different SU{N) groups is identical 
up to a multiple of the identity operator, so it will lead to the same diagrams for the eigenstates. 
Differences follow from the fact that the point where states with higher winding numbers are no 
longer independent depends on N . Considering again the example discussed above of a three-line 
one-loop diagram for SU{2) we find the states, by diagonalising the electric operator, 

(t>i = (123) + (213) + (132) + (321) + (231) + (312) , 

(t>2 = (123) - (213) - (132) - (321) + (231) + (312) , 

03 = (123) -(312) , 

04 = (123) -(231) , 

05 = (213) -(321) , 

06 = (213) - (132) . (85) 
In the one-loop case 05 and 06 vanish, and 03 and 04 are identical. Explicitly, we find 

01 = eClQTrf/) , (86) 

02 = , (87) 

03 = 04 = 3C!QTrC/) , (88) 
05 = 06 = . (89) 

Here are Gegenbauer polynomials ||3^. The state labelled 0i is actually the third eigenstate of 
the Casimir operator of SU{2). In general the nth eigenstate for SU{2) follows from the symmetric 
combination of the permutations of n elements, as other combinations, which contain traces of 
higher powers of U, collapse to lower eigenvalue states. In the example above 03 reproduces the first 
eigenstate. 

In this way we can find eigenstates for fully linked traces (i.e., where the paths have a link in 
common); we can obviously use these results to get eigenstates for unlinked traces by multiplying 
eigenstates of the separate segments. The eigenvalue corresponding to such eigenstates is the sum 
of the eigenvalues of each segment times their length. We also have not yet considered expressions 
containing both the matrices U and . In that case the collection also contains elements where the 
indices of C/'s and C/^'s are contracted with each other, which leads to a cancellation. However, this 
should be a straightforward extension. 



5 The one-plaquette problem 

The above approach is designed for application with a many-body method such as the coupled-cluster 
method [5^. Before doing so we wish to consider the simplest limit, where all plaquettes occur 



independently |36, 37|. Whether these plaquettes contain one, two, or four X variables, Eq. (^), 
the plaquette variable is the only relevant combination. Other variables are trivially integrated over, 
which in principle requires knowledge of the group-invariant Haar measure. For the purposes of the 
present paper knowing the existence of such a measure is sufficient, and all but one of the X variables 
around the plaquette are integrated over. This last remaining matrix we call U. The wave function 
depends on the trace of powers of the plaquette variable only, which can be re-expressed in terms 
of class label, a generalisation of the group characters x{U). It is slightly simpler to investigate the 
U {N) problem first; the SU {N) case can be simply embedded by imposing the constraint det(U) = 1. 
We wish to analyse the Hamiltonian in various ways, both in terms of the trace variables 

ep = -Tv([/p), 
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(we divide by p for future convenience) and in term of the eigenvalues of the matrix [/, which are N 
unimodular numbers of the form e*'^'. The latter parameters are often used in a group-theoretical 



context, see in particular Weyl's seminal work |35]. The most complicated part of the problem arises 



in the evaluation of the action of the electric operators, or, in group theoretical terms, the quadratic 
Casimir operator, on a function of either the ^'s or ijj^s. The potential is rather trivial being the sum 
of traces of U and U'^ . 

If we start with a function F of the ^'s, we find that the action of the electric part of the 
Hamiltonian can be expressed in terms of these trace variables as 

EfFm) = ^^(^TrlC/^+^l'jF,, 

p,q ^ ^ 

/ Ar2 1 \ 

P \ 9=1 / 



(90) 



which can also be generalised to multi-plaquette problems. (Here the notation denotes the 
derivative with respect to ^p, and multiple indices denote multiple derivates.) This relation follows 
from an identity for the generators of U{N), 

E '^ij'^kl = "^^il^jk ■ (91) 
a=l 



The first part of the right-hand side of the expression ( pjOj ) is in principle a complicated function of 
the .^'s if p + g > A^, and we see that the action of the Hamiltonian will be highly non-trivial in these 
variables. The same holds for the trace of the inverse of U contained in the potential. 

Thus solving the one-plaquette problem in this way is difficult, since it involves re-expressing 
Tr[f/P] where p < or p > A'^ in the original trace variables. We shall show that it is easier to use the 
eigenvalue variables, and to perform at the same time a simple transformation which maps a problem 
which is symmetric under the interchange of the eigenvalues (i.e., a bosonic problem) onto a fermionic 
problem, as discussed in the next section. This is similar to, albeit slightly more complicated than, 
what happens for Yang-Mills theories on a circle [^]. 

5.1 Angular representation 

The trace of a matrix is invariant under unitary transformations, Tr[?7"] = Tr[(yf7y~^)"]. Therefore 
we can diagonalise the unitary matrix U in the one-plaquette problem. Prom the N'^ degrees of 
freedom for U{N) only N are relevant to us, 

U = diag (e^'^i , e'^f"^ , . . . , e'^f"^-' , e'^^ ) , (92) 

where for SU{N), ^ipi = 0, since detU = 1. The eigenvalues exp{i'ipi} of the unitary matrix U lie 
on the complex unit circle, as easily follows from the properties of Hermitian matrices by the Cayley 
transformation. 

We use Eq. (|90|), and expand the trace variables in terms of eigenvalues. 



1 ^ 
p ^ 



and assume /({^}) = g{\^Y). It now becomes clear why we have absorbed the factor Xjp in the 
definition of the trace variables, the Jacobian matrix is of the Vandermonde form 

^^j^ = ie^P^^ . (94) 
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This allows us to re-express the action of the electric energy on a general function g{{ilj}), 

E'g={{J-Xg,,),.{J~')jp (95) 

After some formal manipulations involving Vandermonde matrices, which are discussed in the ap- 
pendix, we find the simple result 



2^-2 

i 



(96) 



a Laplace-Beltrami operator in a curvilinear coordinate system. The quantity S is proportional to a 
slightly generalised Vandermonde determinant, 



f^_-^N{N-l)/2 
^_-^N(N^l)/2 



^i(N-l)-4^i/2 gi{Ar-3)»/'i/2 
gi(Ar-l)V2/2 gi{Ar-3)»/'2/2 



gi(Af-l)^jv/2 gi(Af-3)^jv/2 
i i<j 



^-i(N-3)-4^i/2 g-j(Af-l)Vi/2 
^-i(N-3)iP2/2 g-j(Af-l)V2/2 



g-j(Ar-3)i/<jv/2 g-i(Af-l)V'jv/2 



2^^n^i"((^*-^^-)/2) . 

i<j 



(97) 



The intermediate expression masks the fact that H is a real quantity, which is more obvious from the 
first and last expressions. In related group-theoretical literature, especially those following Weyl's 



original analysis |38| one often meets the quantity 



A 



J{N~l),Pi gi(iV-2)Vi gi(Af-3)i/)i 
J(N~l)tP2 J(N-2)il)2 pi(Af-3)i/'2 



gi{Ar-l)^iV gi(Af-2)Viv ^i{N-3)^N 



i<j 



(98) 



rather than H. These are simply related by 



N(N~l)/2 



exp -i(iV-l)/2j^Vi A 



(99) 



The Jacobian J = r? vanishes if two eigenvalues Xi = expitpi are identical. The Jacobian is 
antisymmetric under the permutation of eigenvalues, which corresponds to the reflection symmetry 
on the manifold. For U{2) we can analyse these properties in full. A unitary matrix U £ U{2) can be 
represented as U = e**(n'^ -|- in • o"), with (no, n) real and ng -|- = 1. The real phase \I' is arbitrary, 
but is given explicitly as ^' = —iln[detU]/N . A simple diagonalisation shows that this has the two 
eigenvalues x± = e^'^{nP it i\n\), which are identical for points on both poles, n = (±1,0,0,0). At 
these points Tr[C/] and Tr[C/^] are invariant under first-order changes in the eigenvalues. This would 
correspond to adding a small 6n part to n, but the change in uq is second order in this variable, 
and the imaginary parts cancel to first order. Thus near the poles one of the two coordinates is 
superfluous. 

This becomes even clearer once we realise that the trace variables are equivalent to the parametri- 
sation through the use of diagonal matrices, with the two parameters nP = [e^^'^ + e*'^^)/2 = 
e** cos(^i — V'2)/2 and it' = i(e*^i — e^^'^)/2 = e*"^ sin(^i — V'2)/2. Clearly all points with ipi 
equal to V2 correspond to identical eigenvalues, and for those points one coordinate is superfluous. 
That the behaviour is like that of polar coordinates in three dimensions is not obvious from this 
discussion, and can only be found in an explicit calculation. 
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We have not yet discussed the integration measure, but it can be shown (see Ref. and the 
discussion in Appendix 0) that it takes the form we would expect from the kinetic energy operator, 

(^ii ^2) = l^^imrMm)^'U^4^^- (^^o) 

Here we have not ordered the eigenvalues of U in any way. Since an interchange of the eigenvalues 
corresponds to the same class of matrices, the wave function must be invariant under this interchange. 
Actually, it must be symmetric, and the wave function is that of a system of bosons on a torus. 
However, the complications of norm and kinetic energy suggest that it would be beneficial to make 
the change of functions 

f (W) = E^iii;}), (101) 

to eliminate the first derivatives from the kinetic energy, and the integration weight from the norm. 
As usual this leads to interesting boundary conditions on the function, since they must vanish where 
S vanishes, and since H vanishes linearly and is antisymmetric under interchange of any pair of tp^s, 
the wave function $ is antisymmetric in the same way. This transformation thus turns the bosonic 
problem with a complicated integration measure into a simpler fermionic problem. 



a=l 




N{N^ - 1) 
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Even though it would appear to be quite complicated to determine the eigenfunctions of the quadratic 
Casimir operator (i.e., the square of -E"), we show in Appendix |A| that we can label each eigenfunction 
with a partition [/] = [/i, . . . , /at] 



,i(/i+(Ar-l)/2)^i „i(/2 + (Af-3)/2)Vi ... pi(/jv-(Ar-l)/2)^i 



,i(/i + (Af-l)/2)Viv pi(/2 + (Af-3)/2)Viv ... pi{/iv-(Af-l)/2)^iv 



(103) 



where [/i, . . . , Jn] is the standard U{N) representation labelled by a Young tableau, satisfying fi > 
fi+i- Note that we do not require that all /'s are positive or zero. More details on the representation 
of the symmetric wave function as a one-dimensional fermion problem can be found in Appendix ^ 
where we prove that all excitations <I>[j] can be represented in terms of trace variables. 

Even though this is all quite elegant, we are really only interested in SU{N), rather than U{N). 
Fortunately it is very simple to show that 

^^^^=e'^,f.'^q,^^^(^{^^_^}) , (104) 

where 



N 
J 

is the coordinate describing the phase of the determinant of U. If we interpret \I' as a centre-of-mass 
coordinate, this shows that the eigenfunctions of the kinetic energy (electric energy) separate into a 
trivial centre-of-mass part, and a part that does not depend on this coordinate at all, since it only 
depends on the distance between the centre of mass and the individual particles. This can also be 
written, as we have shown explicitly in the case of H, in terms of the relative angles {ipi — ipj)- This 
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is clearly a simple consequence of the separability of the kinetic energy in centre-of-mass and relative 
coordinates. 

The only well-known problem arising from this separability is that not all functions <I>jj] give rise 
to different functions of the relative coordinates. Thus the solutions where [/] is a partition of 
length — 1 form a complete set of polynomial solutions, as a partition of length can be mapped 
onto partition of shorter length since det U = 1, as we can see from the relation, 

$[/,+n,.../,+n]=e*"^*1>[/„...J,] , (106) 

in which n S Z. 

This knowledge is useful when constructing the full spectrum and all the eigenvalues of the electric 
part of the Hamiltonian, 

N{N'-1) \ _ , 
Dn ^ I ^[/] - %^[/] , (107) 



where the eigenvalues -E^^j follow from Eq. (102), 



N{N^-1) ^ (/, + (iV - 2fc + l)/2)^ , 1^^^ 1 (^\' 



[/I - 24 2 ^8^---^ 

fc=i \fc=i / 

N-^^-^ 1 ^"^ 1 

i=l k=l k<l 

Even though one can gain considerable insight from this separation, we shall generally prefer to work 
with systems described in terms of the "single-particle" coordinates Vi> rather than the centre-of- 
mass and relative coordinates, since it is very complicated to deal with anti-symmetry in these latter 
coordinates. This is well known in the translationally-invariant treatment of the finite many-body 
problem, as can be found in standard nuclear physics text books (see, e.g., Ref. [^). 

As seen above there are no problems when dealing with the electric energy. It is very naturally 
expressed in terms of the momentum operators pi = — ^ as He = \ YliPf- The picture changes 
once we look at the magnetic term in the Hamiltonian which is easily expressed in terms of the 
eigenvalues of U ^ the variables conjugate to p,, as 

Hm = 2N -Tv[U]-Ty[U-'^] 

N 



2^(1- cos ^i) . (109) 



i=l 



There seems to be no easy way to separate this part of the Hamiltonian in terms of a relative 
and centre-of-mass part, and further investigation shows that the commutator of P = ^jPi, the 
centre-of-mass momentum, with the magnetic term is equally non-trivial. 



P, ^ cos tpi 



i^sinVi , (110) 



and thus the potential energy term also breaks translational invariance. The Hamiltonian does not 
commute with the constraint P = 0, which we need to impose on the SU{N) problem, since ^ is an 
unphysical variable. If we deal with this through the standard technology of Dirac quantisation Q, 
where we first classically replace the Poisson bracket by the Dirac bracket 

{A,B]y,b = {A,B]-{A,P]{^,B] + {A,^]{P,B] , (111) 
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Figure 4: An example of the multiplication of the state [2, 1^] in SU{5) with Tr[U] and Tr[[/^^]. This 
corresponds to the addition or subtraction of a square, such that a higher row is always larger than 
or equal to a lower one. A full column can be added or subtracted without changing the state. This 
is used to obtain a canonical form with the last column empty. 



defined so that the Dirac brackets of the constraints ^' = and P = with the Hamiltonian are zero, 
and then quantise the theory by replacing the Dirac bracket by a commutator, we find the relations 
(we use a tilde to denote the variables obeying the Dirac bracket) 

[Pi,Pj] = , 

\pi,i^,] = ■ (112) 

These relations can be satisfied by ordinary commutators if we relate the SU{N) variables ■0 to the 
variables for U (N) by -0 = ^0 — equating pi to pi. In that case the realisation of the Hamiltonian 

i7 = 2^p2^2A^[l-cos(0,-^)] , (113) 

i 

is explicitly translationally invariant, and we can separate all the eigenstates into a ?7(1) part and an 
SU{N) part. An additional factor 4 in the kinetic energy arise from the four separate links around 
a plaquette. The Hamiltonian is scaled with g'^a^^, such that it maps on an A^-body problem with 
unit masses and we introduce the coupling constant A = g^^, for D = 3. 

5.2 Tr[f/] as a raising operator 

We can of course work directly in terms of the angular variables as discussed above. For some cases 
it is actually useful to consider Tr[C/] and its complex conjugate as raising and lowering operators, 
which act in a simple way on the Young tableau labelling the symmetry of the wave function, 

N-l 

i=Q 
N-l 

Tr[f/-']^[/] = ^^[fu-J.-i,-jM-uM . (114) 

where /at = can be imposed, as discussed above. If at any stage in the calculation we end up with 
In 7^ 0, we can use Eq. ( |106D to make it zero. In summary, we need to use the three conditions 

(/l,/2,--- ,/7V-l,-l) = (/l + l,/2 + l,--- ,/^-l + l,0) , 

(n + /i,n + /2,--- ,n + /Ar^i,n) = (/i, /2, • • • , /at-i, 0) , 

(/i,/2,--- Ji,--- Ji,--- Jn-i) = . (115) 
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Even though this interpretation of the traces is quite appeahng, it is not totally straightforward to 
use, since there is no cyclic vector; when acting on the state labelled by 



[/] = [0^-i] , (116) 

both terms lead to non- vanishing results. This is intimately related to the existence of both covariant 
and contravariant representations of SU{N) occurring at the same kinetic energy, which are mixed 
by the magnetic terms. An example is shown in Fig. |^. 

5.3 The harmonic approximation 

In the weak-coupling limit we can make the harmonic approximation, where we approximate the 
magnetic potential by its quadratic expansion, and we can ignore the periodicity of the ip variables 
since they oscillate near their equilibrium values. The resulting potential, 

N 



Vim) = x^{i^^-^? 



i=l 



(117) 



is modified by confining the centre-of-mass in the harmonic oscillator (HO) potential 1^(^') = ^N"i!^. 
We have reconstructed here the well-known separability of the many-body harmonic oscillator Hamil- 
tonian into a centre-of-mass (CM) and relative part, where the CM part is in a harmonic oscillator 
state with frequency N times larger (since the mass for the CM mode is also N times the single- 
particle mass) than the single particle one (and see, e.g., Ref. |pT|). 

Amongst the spectrum of the A^-fermion HO problem we have states with the CM in the m = 
0,1,2,... states, and we thus reproduce each relative state with all of the CM states. If we can 
determine all those states with the CM in its ground state, we immediately know the degeneracy of 
the SU (N) spectrum. So we need to identify those states from amongst the degenerate multiplet 
with energy {n + N/2)huj. Let Un be the degeneracy of this (n-th excited) energy level, which is given 
by the number of distinct partitions of n into N positive numbers: 



z/---' =Pjv(n) = Yl 1 • (1^^) 



Ai>A2>- >Aiv>0;E Ai = 



We use the fact that harmonic oscillator wave functions are exponentials times a polynomial. Clearly 
the centre-of-mass polynomial can only arise from a linear combination of the polynomial parts of 
the various degenerate wave functions. For a CM ground state this polynomial part is a constant, 
and thus lies in the null space of the centre-of-mass momentum operator 

N „ 

p = yi^ . (119) 

This operator maps the n-th degree polynomials onto the (n — l)-th degree polynomials. The di- 
mension of the null space is the difference between the dimensions of the domain and range of the 
operator P. Therefore the degeneracy of the n-th relative motion, i.e., SU{N), eigenstate is given by 

i/f^(^) = P^(n)-P^(n-l) . (120) 

For A^ = 2 

„ , X r?^ + ii /X 

P2{n)= , (121) 
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where [ • ] denotes the integer part of a number, and thus f„ = (^n,odd) which corresponds to a 
set of non-degenerate, odd wave functions for SU{2). For SU{3), S'C/(4), and SU{5) we find the 
degeneracies (the first state always occurs at i = N{N — l)/2, since ah lower energy functions are 
not antisymmetric). 
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(122) 

Notice the absence of the first excited state; i^^l^iy2+i ~ ^' since the only state at that position is 
a centre-of-mass excitation built on the ground state. Note that our result differs from the result in 
Ref. [^2|, where the fact that only even functions in the radial coordinate are allowed was ignored. 
The SU{3) spectrum is 



^2V2\{3k + 2n) , (123) 
5 9 

for n = 0, 1, 2, 3, . . . and k = 1, 2, 3, ... . 
5.4 Large-A^ limit 

The large- limit of QCD is one of the intriguing open problems in the theory, and one might wonder 
whether our result can shed some light on the much less taxing problem of the large- limit of the 
one-plaquette problem. In that case we have to study the Hamiltonian 

N 

i/ = 2^pf + 2Aj;(l-cos(^,-vI/)) . (124) 

i i=l 

The leading term is simply obtained from the Hartree-Fock approximation, with Hartree Hamiltonian 

n = 2p^ + 2A(1 - cos V) , (125) 

and we can (if we wish) evaluate the Hartree-Fock energy for this state using the solutions of the Math- 
ieu equation. The boundary conditions are strict periodicity, since the all the circles parametrised 
by V'j are connected, and end up at the same physical point. 



Such results are well-known from many-body theory |43|. In principle the next order in 1/A^ would 
follow from an RPA calculation. However, there are some problems. The interaction is explicitly 
A^-dependent, since ^ depends on A^, and the residual interaction is rather cumbersome to deal with. 

5.5 Explicit solutions 

Let us now investigate the spectra of a few of the relevant gauge theories. The method we use to solve 
the problem is to work in a basis of eigenstates of the electric Hamiltonian, and evaluate the action 
of Tr [/ and Tr C/^ on these states. As can be seen from the example for SU (3) shown in Fig. ^, there 
are several different coupling mechanisms. Nonetheless, this method is straightforward to implement, 
and leads to fully converged results up to large values of the effective coupling constant l/g^- 
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Figure 5: Examples of the action of Tr U (solid lines with arrows) and Tr[J7^] (dashed lines with 
arrows) on the states of SU{3), labelled by the partition [/i, /2]. 




Figure 6: The lower end of the fully converged U{1) one-plaquette spectrum. The lines starting 
from 4n^ are the periodic solutions, and the lines starting from (2n + 1)^ are anti-periodic ones. The 
hashed areas denote the existence of a continuum of Floquet states at these energies. 
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5.5.1 The U(l) one-plaquette problem 

First we look at U{1). This is something of a special case, as it does not fit into our general SU{N) 
framework. However, the relation is straightforward. With the parametrisation U = e**^, the electric 
operator is id/dcf). If we use the maximal tree formulation, we find that the action of the electric 
energy gets weighted by a factor of four, essentially one for each side of the plaquette. The magnetic 
part is [2 — exp{i(j)) — exp(—i(f))] = 2(1 — coscp), and hence the Schrodinger equation becomes 



2a 



d2 



cos ( 



(126) 



Upon substitution of 2z 



a = 2{Ea-l)/g^, q 



-2/g'^, we find the Mathieu equation, 
a — 2q cos2z)y{z) = . 



(127) 



However, we still need to analyse the boundary conditions on ^. It is well known from Floquet 



theory [44| that solutions of a differential equation with periodic coefficients are quasi-periodic. 



(128) 



where ^{(p) = ^{(j) + TT), and v is a, in principle complex, exponent. The only normalisable solutions 
are those with a real exponent, and those are found bracketed between the (periodic) cer{z,q) and 
the ser+i{z, q) solutions [^], as in Fig. These results are analogous to the band structure in metals, 
due to the periodic potential in which the electrons move. We define a scaled energy e as 



^ N 
E = ^e + — 
la g'^a 



(129) 



which we use for representing the numerical results. 

For the period states we can attack the problem numerically by expanding the Hamiltonian in 
terms of the eigenstates, exp(m(/)), of the electric operator. 



El 

n=0 
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n=l 



n) 



n) 



4=^ + il(»l 
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2a 



g'^a 



2g'^a 
1 



{|n)(n + l| + 
{\n){n + l\ + 



n 



\n ■ 



+ l)(l + 5o„)(n|} (130) 
-l)(n|} (131) 



where the sum starts at zero for ce2r solutions and at one for se2r solutions. This matrix equation 
serves as a relatively efficient way of determining the characteristic values a and b of the Mathieu 
equation, via straightforward numerical methods. 



5.5.2 The SU (2) one-plaquette problem 

The SU{2) matrices have two complex conjugate eigenvalues, and can be parametrised as Tr(C/) = 
2 cos (f), where cf) = (t/'i — '^2)) and see Eq. (p2[). If we absorb the factor H from Eq. ( |97| ) into the wave 
function (H = 2 sin(i?!i/2)), 



The Schrodinger equation, which takes the form 



r 

2a 



a2 



cos(0/2)) 



(132) 



(133) 
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Figure 7: The lower end of the fully converged SU (2) one-plaquette spectrum 

is again the Mathieu equation. In this case the interpretation is different, even though we still have 
2z = (p, but (p plays a different role. As a function of the periodic quantity TV [/, $ must be strictly 
periodic. Due to the anti-symmetry induced by the Pauli principle, the only non-singular wave 
functions ^ are the odd Mathieu functions sei{z,q),se2{z,q),- ■ ■ , for which a = br{q), where the 
index is allowed to take odd values, since the 27r periodicity on [— vr, vr] is guaranteed by the solutions 
being odd. For q = the solutions reduce to = sin^^c/), and since H oc sin^cj), ^'(— tt) = ^'(vr). 
The solutions are shown in Fig. 0. 

In terms of the eigenstates |n) of the electric operator the Hamiltonian reads 

H = y\n) f '^^'' ^ + 4-] {n\ - -^{\n){n + l\ + \n + l){n\} . (134) 

We see that the differences between the U{\) case and the SU{2) case lie in the diagonal part of the 
Hamiltonian, respectively and (n -|- 1)^ — 1, which result from the electric operator. 

5.5.3 SU{3) one-plaquette problem 

The SU{3) one-plaquette problem has two independent angular degrees of freedom. From the U{3) 
parametrisation we see that these can be chosen as (pi = ^{21^1 —ip2 — V's)] (^2 = ^("V'l + ^^^2 — V'a)- 
The remaining combination, ^{—ipi — ^2 + 2V'3) = — {(pi + '^2), is linearly dependent. 
The Hamiltonian in these angular coordinates becomes, using Eqs. ( |102| ) and ( |109|) , 



9^ 
2a 



3\ d(Pl d<pl^ d(Pid^2 ]^ d 



2-1^ + - 3 r + 74(3 - COS</)l - COS</)2 - COS(01 + <P2)) 



^ = E^ , (135) 



where the boundary conditions are determined from the anti-symmetry of ^. The wave function 
vanishes if two eigenvalues of the matrix U coincide, (pi ^ <p2 ~{4'i ~^ 4'2) 7^ <Pi- All variables 
are defined modulo 2tt, as they correspond to identical eigenvalues. We can restrict ourselves to the 
domain —(pi — (p2 ^ (pi , (pi ^ (p2, and (p2 < 2tt — (pi — (p2- However, many other choices of boundary 
conditions are possible, which follow from either permuting the variables, or translating the domain 
by 27T. The (pi Cg) (p2 plane is divided into copies of the fundamental domain by the lines (pi — (p2 = ^itt, 
(pi + ^(p2 = n2TT, and ^(pi + (p2 = n^Tr, where G • • • , —2, —1, 0, 1, 2, • • • , as shown in Fig. ^. 
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71 271 371 471 

Figure 8: The domains for and 02, in the antisymmetric wavefunction 
The potential can be simpUfied to 

1 , A I A , (A , J. \ A (Al + 4>2 

1 + COS 01 + cos 02 + COS(0l + 02 j =4 COS — COS — COS . 

The exphcit angular Schrodinger equation is useful for the harmonic approximation in the weak- 
coupling limit, however, for the strong-coupling limit and the intermediate region in coupling con- 
stants, the Hamiltonian in terms of the electric-operator eigenstates /2) is more useful 

' (/l,/2 

\ nil I 

o<h<h 
1 



■{|/l,/2)(/l,/2 + l| + |/l,/2)(/l + l,/2| 



+ l/l + l,/2)(/l,/2| + |/l,/2 + l)(/l,/2| 

+ |/l,/2)(/l + l,/2 + l| + |/l + l,/2 + l)(/l,/2|} , (137) 

where all the states |/i, /2) where the condition < /2 < /i does not hold, are identical to zero. This 
problem does not seem to have a closed-form analytical solution, and we have to resort to numerical 
methods to obtain the spectrum for arbitrary values of the coupling constant g, as in Fig. |9|. 

5.5.4 Higher orders 

We conclude this section by showing results for the gauge groups S'C/(4) (Fig. [l^) and SU{5) (Fig. |ll]). 
The most interesting aspect of these results is the presence of (what appear to be) real level crossings, 
which could be a reflection of some hidden symmetry in the problem. It was verified that the distance 
at these crossings was equal to zero within the numerical accuracy. 

The region of the coupling constant shown in Figs. |^-|ll| is insufficient to see convergence to the 



harmonic approximation (122). We have checked numerically, for g > 50, that these results are 
correct. 

6 Wave functionals 

The results for the one-plaquette problem have more consequences for variational wave functionals 
than one might suspect. If the trace of the one-plaquette matrix is used, the wave functional is a 
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Figure 11: The lower end of the fully converged SU{5) one-plaquette spectrum 

function of the group characters only. The wave functional, which is the sum of of one-plaquette 
functions, 



{{^l^a^m= Yl ^ [{^c,^}l-,') , (138) 
plaqucttes a 

naturally leads to the sum of one-plaquette problems, leading to the energies which are the sum 
one-plaquette energies. However, the product wave functional 

{{^a^m)= n F{{^a^}l-^') , (139) 
plaquettes a 

also leads to the same result. Note that the exponentiated potential, exp{HiHM}, where k is a 
parameter to be determined by a variational calculation, as a simple ansatz for the correlated wave 



functional |4(:] is a particular choice in this class of wave functionals. The absence of correlations 
between nearest-neighbour plaquettes, generally present, follows from the symmetry in the angular 
variables, 

^^(...^^...^....) = ^^(...^....^....) , (140) 

where 

F(Vi---Viv) = X]cA^A(^i---V'iv) . (141) 

A 

Therefore the product term from the electric operator vanishes, 

DN^^'xii^al ■ ■ ■ V'aAr)^A'(V'/31 " " " V'/3Af) = ^A' (^/31 ' " " Pn)D N^'xii^al ' ' ' '^cin) 

+ ^'x{'^al---i'aN)DN^'x'{i^pi---i'l3N) , (142) 

where a and (5 are plaquettes containing the link /, and = J^\. The differential operator Di^ 
contains both sets of angular operators 5^. = 9^^. + d^^.. Therefore the Hamiltonian, acting on the 
product wave functional |n) also reduces to the sum of one-plaquette Hamiltonians. 
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7 Conclusions 



With the understanding of both a gauge fixed Hamiltonian, and the one-plaquette problem we would 
now like to use this in a method designed to determine the ground-state wave function and excitations 
from the Hamiltonian. Since there are, in principle, an infinite number of degrees of freedom, some 
approximation is required. However, we want to do this in a systematic manner, where we can improve 
the order of approximation as required. The coupled cluster method (CCM) |^, a many-body 
technique, is suitable for such calculations. It has previously been applied lattice field theories, [p7| , 
and in various guises to gauge theory although in a different context. It has also been applied 



to various spin-lattice problems in magnetism |4£]. 



Central to the CCM is the parametrisation of the ket-state wave function as an exponential of 
excitation operators 

|*)=e^|^') . (143) 

The bra-state contains the inverse exponential such that the Hamiltonian functional only contains 
linked terms of the excitation operators, due to the implicit similarity transform He^ . The 
projection on the model state |<1>) is due to the annihilation operators S, 

= (<I>|(l + 5)e-^ , 
(^1 = (^le-^e-^ , (144) 

where the first form is used in the normal coupled cluster method (NCCM), and the second form 
in the extended coupled cluster method (ECCM). The normal coupled cluster method does not 
reproduce the correct scaling limit when the excitation operators are obtained from the strong- 
coupling expansion. 

However, this and the extreme difficulty of implementing the ECCM are related to the functional 
implementation of the CCM, where an operator is simply chosen to be the multiplication of the wave 
function with a function. We will only consider the operatorial CCM, where the creation operators 
excite from the model state to any arbitrary state, and do not act between different excited states 
on overlapping lattice-site configurations 

S = Y,<^i\I)m ■ (145) 
I 

This allows us to implement the ECCM, which is better at describing a system at large global changes 
away from the model state 

The central question is the form of the states |/) in which to expand the Hilbert space. This is 
the question discussed in the first part of this paper. For the colourless sector we have to use closed 
contours, which are traces over products of X variables, since only these variables are invariant under 
gauge transformations generated by Gauss' law. The traces over X variables are automatically closed 
contours as each contour associated with an X variable starts and ends at the origin. 



Our Hamiltonian is different from the naive Kogut-Susskind Hamiltonian used in Refs. |31,32,5C], 
and this should have an important effect on the results we hope to obtain in our future work. The 
price we have to pay is that our Hamiltonian has a preferred direction, and we have lost explicit 
translational invariance. It is probably worthwhile to investigate different, more symmetrical, choices 
of the maximal tree. It would seem that the current choice should be very good to consider the 
interaction between fixed sources, which clearly break translational symmetry, but it is less obviously 
optimal for a study of the vacuum. 
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A One-dimensional many-fermion systems and their representa- 
tions 

As the SU (N) one-plaquette problem can be mapped onto a one-dimensional (A'' — l)-fermion prob- 
lem, we analyse the fermion problem in some depth to determine the spectrum and the action of a 
potential in coordinate representation. Usually fermion problems are dealt with in operator form, 
here we use the knowledge of (anti) symmetric functions to formulate the problem. 

An n-fermion wave function is in terms of single-particle orbitals 4>i{x) given by the fully anti- 
symmetrised Slater determinant 



(Xi,... ,Xiv|*) 



(pi{xi) 



4>2ix2) 



(t>3{xi) 



(pliXn) hiXn) hiXn) 



(f>n{xi) 
(I>n{x2) 

(f>n{Xn) 



(146) 



If 4>m{x) = w{x)Pm{x) are chosen to be the lowest set of single-particle orbitals, in which w is the 
weight function, which is the square root of the integration measure, and Pm is the m-th orthogonal 
polynomial given that weight, then the wave function ^ is the ground state wave function $o- Only 
the linearly independent parts of each column in the Slater determinant contribute, and these are just 
the highest powers in each polynomial Pm{x) x"^. Therefore the ground-state Slater determinant 
reduces to the Vandermonde determinant 



(x|$o) 



1 



JJ'u;(xi) -xi) 

1=1 k<l 



(147) 



As examples we deal with the infinite square well 
one-particle wave functions are given, respectively, by 



-7r/2,7r/2] and the harmonic oscillator. The 



{'?!'}s.w. 

{(t>Wo. 



—= sin 2nx; 



—= cos(2n -I- l)x 

"I'^Hnix) . 



(148) 
(149) 



The sines and cosines can be re-expressed as polynomials in sinx with the weight cosx. The infinite 
domain of the harmonic oscillator is of no consequence as the weight damps the wave function at 
large distances. The non-interacting n-fermion ground state consists of fermions in the lowest n 
levels. We find 



(a;|$o)s 



(a;|$o)h.o. 



21-/2] [n/2]!" n 
' I I cosxj I I (smxfe — smx;) 




k<l 



-xl/2 



H^iXk 



(150) 



(151) 



k<l 



where the normalisation of the square-well wave function follows from De Moivre's formula, and the 
normalisation of the harmonic oscillator wave function is a direct consequence of the normalisation 
of the Hermite polynomials Hn- 

We now consider an arbitrary excited state ^x. Since the anti-symmetry is present by the ground- 
state wave function, the wave function of an excited state can be factorised as a symmetric function of 
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Figure 12: The labelling of excitations in a fermionic system as non-distinct partitions, where any 
partition A1A2 • • • for which A, > Xj if i < j represents a unique excited state \n{n — l)(n — 2) • • • ) — > 
|(n + Ai)(n — 1 + A2)(n — 2 + A3) • • • ), where the labelling of states is the distinct partition. 



all the variables times the ground-state wave function. For any finite basis, this symmetric function 
must be a symmetric polynomial 



{x\<i>x) = Sxix){x\<l>o) 



(152) 



Now the possible excitations follow from the theory of symmetric polynomials, and can be labelled 
as a distinct partition A = (Ai, A2 • • • , A„), where Aj is the label of an occupied level, in ascending 
order. 

The symmetric polynomials are usually labelled independently of the degree of filling n by the 
non-distinct partition f = X — 5 where 6 = (n,n — 1, • • • ,2, 1) such that fi < fj for i > j. The 
non-distinct partition / yields the excitation with respect to the Fermi level (see Fig. |l^). 

Generally, the symmetric functions S\ are not of a definite degree and they have a normalisation 
that depends on n. So A refers to the leading order polynomial S\ = cs\ + • • • . These polynomials 
sx are the Schur functions |51|, 



i 

52 = ^Xi+^XiXj , 

i i<j 

Sl l — ^ ^ XiXj 
i<j 

53 — ^ ^ ~l~ ^ ^ {Xj^ Xj -\- XjXi) + ^ ^ XiXjXk , 

i i<j i<j<k 

S2,l = '^{xfxj + x'jXi) + 2 ^ XiXjXk , 
i<j i<j<k 

Sl,l,l = ^ XiXjXk ■ 
i<j<k 



(153) 



For SU{N) one-plaquette problem the excited states are symmetric polynomials in the eigenvalues 
of the matrix U £ SU{N). These polynomials are the Schur functions, S\ = s\, and, as the 
eigenvalues lie on the complex unit circle, the symmetric polynomials are properly normalised. The 
eigenvalues of the electric operator follow directly from the definition. The two remaining problems 
are the action of the magnetic term on these states and how, and if, to express the symmetric 
polynomials, in the eigenvalues, again in terms of trace variables. For this purpose we introduce 
three representations: the Schur functions s, the monomial symmetric functions m, and the so-called 
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trace functions t, 



l+Ajv-i+i 



det X- 
sx = ^— T 



det X* 



Perm{jj) 

ta = nl^E^j • (154) 

The numbers A denote a partition, and the a label can be related to a partition by 

A = •••2"n"i) . (155) 

For instance, a = (1020 • • • 0) corresponds to A = (331). 

Prom the relations between the different representations one can show that every excited state 
can be expressed as a product over traces. An excited state s\ of the SU{N) one-plaquette problem 
is a partition of length less than N. If we use a lexicographical ordering >~ |5l|, it means that S\ 
contains only monomials where k ^ X. It is easy to see that a trace function only contains 
monomials where n ^ u. Therefore, an excited state ^\ from Eq. (pOD, can be expressed in terms 
of the — 1 independent trace variables, where 

taia^.-a^.i = Tr[[/]-^ Tr[C/2]°2 . . . Tr[C/^- Y^-^ • (156) 

We can thus express the anti-symmetric wave functions in terms of the symmetric trace variables. 
One can easily construct the transition, or Kostka, tables |^l| for transformations between s and m, 
as well as t and m. 

Even though it is relatively straightforward to use trace variables for SU{2), it turns out to be 
rather involved to extend these relations to higher-dimensional gauge theories. For example, with the 
relations above we can determine the integration measure for trace variables. For SU{2) we require 
the relations up to the partitions of 4, while for SU (3) we require similar relations for the partitions 
of 6, 

SU{2) : E'^ = l-tj , (157) 
SU{3) : = 27 - 5t? + Otita - ^4 + - 1% + (158) 

= 27 - 18iit_i + 4(t? + ti^i) - . (159) 

where tp = Tt\U^]. However, the use for this relation for SU{2>) is restricted as the two complex 
trace variables ti and t2 have complicated relations, and are restricted to a non-trivial domain of 
the complex space C^. Note that there are only two independent angular variables for SC/(3), which 
makes the angular variables much more useful than the trace variables. 

The Haar measure //(C/) reduces to the measure of the trace variables if the integrand depends 
only on the trace variables, 

j dfi{U)F{{tx}) = J Edti ■ ■ ■ dtN-iFi{tx}) = I . . . #^_lF({^A({0^})}) , (160) 

where H is a non-polynomial function of the trace variables tj. 
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A.l The magnetic term 



The magnetic terms in coordinate representation are just the multiphcation of the state with si 
and where 11 • • • 1 is 1^~^. We can apply the Littlewood-Richardson rule for the general 

multiplication of two Schur functions, 



E^Ak^/. • (161) 



However, since si and are relatively simple we can derive multiplication rules from first princi- 

ples, taking into account that the partitions cannot be longer than A*" — 1. As symmetric polynomials 
we find {xj = exp{i'>pj}), 

Tt[U] = XI+X2 + X3 + -- - +XN , (162) 

Tr[?7^"'^] = xiX2---XN-i-\ 1- X1X3 • • • xat + X2X3 • • • xat , (163) 

where we use that X1X2 ■ ■ ■ xn = 1. In terms of the partition A of s^j multiplying with Tr[J7] means 
adding one to each in turn Aj ^ Aj + 1, cancelling all the partitions for which Aj > Xj for j < i, due 
to the intrinsic anti-symmetry. Multiplying with Tr[C/^^] corresponds to Aj — > Aj — 1, with the same 
cancellations. 

This relation can be generalised to the multiplication by Tr[C/*'] for any power p. This reduces to 
adding p to each of the elements A, in the partition A. The partitions should be ordered Ai > A2 > 
• • • > Aat-i > 0, and the non-distinct partitions dropped. 



B Equivalent forms of the electric energy 

In this section we prove that the electric operator in angular coordinates {ipi}, 



D 



N 



1 ^ 
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2 ^ di^f + 2 
1=1 ' 



1 / 1 



N 



d 



(164) 



the SU{N) equivalent of Eq. (102), Can be realised in the trace representation {Tr[C/P]} by the 
SU{N) version of the electric operator in Eq. (]90|) . 

We start of by the relation between a function F of the trace variables and one of the angular 
variables, 



F(Tr[C/],iTr[t/2],... ,^Tr[C/^-i: 



/TV N 



.2=1 



2=1 



N 



1 ^ \ 
~^ i=i J 



(165) 



In order to prove the fundamental relation, we need to know the derivatives of the Vandermonde 
determinant with respect to the angular variables 



A. 



A 5V'. 



N 



(166) 



Writing out the terms in Eq. (|102|) we generate several terms 



+ 



II 




III 




(167) 
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where di = Some tedious, but straightforward, calculations lead to 



I = ^ ' f , (168) 



N-l N-l 



II 



i nUP-"^]F,, + Y,pTr[UP]F,, (169) 



yp,q=l p=l 

' N-l N-l 



^ Tv[U^] Ti[U']F,p, + Y PnU']F,p , (170) 

\p,q=l p=l J 

N-l /p-1 -, / 1 \ \ 

2^ lE2^[^'"']T'[^'J+(^-2"lj^'[^'U^'^ ' ^^^^^ 
AT - 1 

IV = _^Tr[C/^]F, , (172) 

p=i 



III 



which, when substituted back in the original equation, yields the SU{N) version of Eq. (pO|): 

\ P J p,q ^ ^ 

+ E ^F^Tr[C/^'] + - ^ TrlC/^'-^lTVit/^] ) . (173) 



q=l 

This finishes the proof. 
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